Born-Oppenheimer approximation in open systems 
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We generalize the standard Born-Oppenheimer approximation to the case of open quantum sys- 
tems. We define the zeroth order Born-Oppenheimer approximation for an open quantum system 
as the regime in which its effective Hamiltonian can be diagonalized with fixed slowly changing 
variables. We then establish validity and invalidity conditions for this approximation for two types 
of dissipations — the spin relaxation and the dissipation of center-of-mass motion. As an example, 
the Born-Oppenheimer approximation for a two-level open system is analyzed. 

PACS numbers: 03.65.Bz, 07.60.Ly 



I. INTRODUCTION 

Schrodinger equation which governs the dynamics of 
quantum system is among the fundamental equations in 
physics. In principle we can predict all features of a quan- 
tum system at zero temperature by solving this equation. 
However, many problems are hard to handle and the 
Schrodinger equation might not be solved analytically. 
Many approximations have been developed in order to 
solve this problem, including the adiabatic approxima- 
tion [l[ and the Born-Oppenheimer(BO) approximation 
i. 

The adiabatic approximation tells us that for a suffi- 
ciently slowly varying Hamiltonian, if the system initially 
is in an instantaneous eigenstate of the Hamiltonian, the 
system will remain in that eigenstate up to a phase fac- 
tor at later times (for the adiabatic condition and recent 
progress, see Ref. 0, 1, 1, @, 0, 1 S EE El El ) • 

When the 

quantum system evolves adiabatically and cyclically, the 
eigenstate of the Hamiltonian will acquire a Berry phase 
[1 31 ] except the dynamical phase. This concept was gen- 
eralized to non-cyclic and non-adiabatic cases, and the 
non-dynamical phase acquired in this situation was called 
the geometric phase EE EE EE EE EE HE H3 • If a 
quantum system has two sets of variables, say a fast and 
a slow set of variables, the BO approximation can be used 
to solve the dynamics, this can be done by solving the 
fast variables with fixed slow variables first, leading to 
an effective Hamiltonian for the slow variables. Apply- 
ing this approximation to a bipartite particle, an exter- 
nal vector potential is introduced for the slowly moving 
particle due to its fast moving partner [22[. Thus the 
total wave-function can be factorized to a product of two 
wave-functions corresponding to the fast and the slow 
variables. The BO approximation has been widely used 
in many fields of physics and chemistry and proved to be 
a fundamental tool in these research [23(. 

These two approximations arc both for closed quan- 
tum systems. However, due to the unavoided couplings 
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of quantum system to its environment, a realistic quan- 
tum system is in general open and the dynamics of such 
a system can be described by the so-call Markovian mas- 
ter equation. The master equation was derived from 
Schrodinger equation by tracing out the environment 
variables in the Markovian limit 1241 (for the solution 
of this equation, see e.g. Ref. [25l. |26|. l27l. l28l. l29l [30j] and 
for the dynamics beyond the Markovian limitation, see 
Ref.[H M, H HE EE HE HE HE El)- The adiabatic 
approximation has been extended to the case of open sys- 
tems by several authors [U |H, EE H EE HE E3 • 0ne 
of these extensions was done by the effective Hamilto- 
nian approach [45l . 1461 E3]. The key idea of the effective 
Hamiltonian approach can be outlined as follows. By 
introducing an ancilla which is the same as the quan- 
tum system, we can map the density matrix of the quan- 
tum system to a wave-function of the composite system 
(system plus ancilla). Then the master equation can be 
mapped to a Schrodinger-like equation. Since the dy- 
namics of open systems is governed by a Schrodinger-like 
equation, some conclusions for the closed system can be 
extended to open systems straightforwardly. In this pa- 
per, we extend the BO approximation to open systems 
by using the effective Hamiltonian approach. Two types 
of dissipation, the spin relaxation and the dissipation in 
the center-of-mass motion(DissCOM), will be considered. 
The validity condition for the BO approximation is also 
presented. 

The structure of this paper is organized as follow. In 
Sec|H]we generalize the BO approximation to open sys- 
tems with spin relaxations. The case of dissipation in 
the center-of-mass motion (DissCOM) is considered in 
Sec lIIII In Sec lIVI we present an example for the case 
of spin relaxation to show the details of the formulism. 
Conclusions and discussions are given in SecfVl 



II. SPIN RELAXATION 

Consider a spin moving in a magnetic field. Taking the 
spin relaxation into account, the evolution of the density 
matrix for such a system can be described by the Marko- 
vian master equation that generally can be written in the 
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following form (h = 1) [24 



-i[H,p] + C P , 



(1) 



where -ff is the Hamiltonian for the system. In our discus- 
sion, this Hamiltonian can be divided into the following 
three parts 



H = H K + V + H S (B(3),S), 



(2) 



where Hk = Hk(x) — is the kinetic energy, p is 
the momentum operator, V — V{x) denotes an exter- 
nal potential and Hs(B(x), S) represents the coupling of 
the spin to an external magnetic field B(x). Note that 
the coupling depends on the coordinate x through the 
inhomogeneous external field B(x). The second term 
on the right hand side of Eq.JI]) denotes the spin re- 
laxation/dissipation, which is assumed here to takes the 
Lindblad form [48|, 



pa + a 



<7 + <T p), 



(3) 



where a + and a~ are Pauli operators. For the case of 
DissCOM, this term is different [Hj]. We shall discuss 
this problem in Sec lIIII 

Following the standard procedure of the effective 
Hamiltonian approach (3ll . |45| . |47| , we can map the den- 
sity matrix into a wave-function by introducing an an- 
cillary system, and obtain a Schrodinger-like equation 

PES si, 

N 

P->\$)= Pmn(t)\E m )\e n ), 
m,n— 1 



(4) 



where p mn (t) — (E m \p(t)\E n ) is the element of the den- 
sity matrix. {|_E„)} and {|e n )} are the time-independent 
bases for the spin and the ancilla, respectively. Ht is the 
so-called non-Hermitian effective Hamiltonian. For the 
considered moving magnetically-driven spin(Eq.([S])) and 
the master equation (TTJ), the effective Hamiltonian has 
the following form, 

Ht = Hk-Hk + V-V a + H s -H§+L 8 , (5) 

where the superscript A denotes operators of the auxil- 
iary system, which are defined by 



|0 / 



(E n \0*\E„ 



(6) 



and Hk = 



H A = 



H A (x A ), V a = V A {x A ) 
Hg(B(x A ),S A ), Ls represents the spin relaxation, its 
form depends on the coupling of the system to its envi- 
ronment and hence it is a function of the jump operators. 
For spin relaxation/dissipation (Eq.Q), we have 



-JT r T 



(7) 



where r + and r~ are Pauli operators for the auxiliary 
system. Notice that there is only one term in Ls de- 
scribes the interaction between the quantum system and 
the ancilla. Now we are in a position to extend the BO 
approximation from closed to open quantum systems. In 
order to simplify the derivation and make the formulation 
more explicit, we define a new six-dimension coordinate 
by the coordinates of the quantum system and the ancilla 
as r=(x, ix A ). With these notations, the effective Hamil- 
tonian can be rewritten as Ht — H\ (r ) + V T (r) + Hg (r ) 
withifj(r) = H K (x)-H A (x A ), V T (r) = V(x)-V A (x A ) 
and Hj(f) = H s (x) - H^(x A ) + L s . H]j(r) is the total 
spin part of the effective Hamiltonian Ht- We denote 
the right eigenfunctions of Hg(r) by |A^) = |A^(r)), 
(n = 1,2,..., A 2 ) and the left eigenfunctions by (A^| = 
(A^(r)| with the corresponding eigenvalues E n = E n (r) 
for a fixed f. Generally speaking, the total spin part Hj 
is not Hermite and the eigenvalues are usually complex. 
Nevertheless the left and right eigenfunctions still satisfy 
the orthonormal relation (A^JA^) = S mn . We should 
note that when we map the density matrix into a wave- 
function of the composite system, degeneracies may be 
introduced. The occurrence of degeneracies is a signature 
of the presence of symmetry in the composite system [50| . 
In the following, we have restricted our discussion on the 
non-degenerate energy levels. 

We can expand the eigenstate of the effective Hamil- 
tonian Ht in terms of {|A^)} as 



JV 2 

|$) = 5> n (r)|A«(r)>. 
fe=i 



(8) 



Substituting Eq.© into the eigenvalue equation of the 
effective Hamiltonian Ht\$) — E\$>), we obtain the fol- 
lowing equation for <!> n (r), 

H T (n)^ n + F T (n)$ n + £ T (n,m)$ m = E$ n , (9) 



where 



1 



y ' 2M 



[V r -iA^n)} 2 +V T (r)+E n (r) 1 (10) 
A>'(n)=i(A^|V ? A«), (11) 
FT{n)= ~2JI £ «An|VrA£><A£|VM*» , (12) 
T (n,m)=-^L (2(A^|V r -Af n )V ? +(A^|V|A«)) .(13) 

Vf in the above equations denotes the derivative with 
respect to the new coordinates r. We can see from these 
equations that |V^A^) = for all m may lead F T and 
T to vanishing, this results in a complete separation 
of the spin and motion variables in the dynamics, which 
is quite different from closed systems. It is easy to find 
one solution that these condition can be met when the 
dissipation is negligible and the external field is homoge- 
neous. Treat F T and O t as perturbations, we can use 
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the standard perturbation theory to solve Eq.©. First 
we rewrite Eq.© as 

(Ho + W)Q = E$ (14) 

where $, Ho and the perturbation W are defined as 





" $1 " 




~H T (1) 











ff T (2) • 









, Wo = 






<f> N 2 _ 







■ H T (N 2 )_ 



This is the validity condition for the zeroth order BO 
approximation. 

In this section, we have extended the BO approxima- 
tion to an open system with spin relaxations. The zeroth 
order approximation was defined as the regime where 
its effective Hamiltonian can be diagonalized with fixed 
slowly changing variables. This generalization is avail- 
able for dissipations in all variables except coordinates 
(or center-of-mass motion), which we shall study in the 
next section. 



W 



F T (1) 
O t (2,1) 



T (1,2) 
F^{2) 



O t (N 2 ,1) T (N 2 ,2) 



O t (1,N 2 ) 
T (2,N 2 )) 

F T (N 2 ) 



(15) 

By making use of the time-independent perturbation the- 
ory, we can obtain high-order approximations for Eq.© 
or (|14j) . The zeroth order approximation solutions 



b 


,$g(2) = 


■ " 

*g(2) 











are given by the eigenvalue equation 



iP(n)<I>L u >) 



(16) 



(17) 



kp in the above equations represents a wave-vector de- 
fined by kp = (k,ik A ). This zeroth order solution in- 
dicates that in the six-dimension space, the system and 
ancilla move in an external potential V(r) as a whole sys- 
tem and evolve according to the free Hamiltonian H T (r). 
The interaction between the spin and the magnetic field 
together with the spin relaxation make the two systems 
feel a complex vector potential which would in turn 
affect the center-of-mass motion of the particle. 

From these zeroth order solutions, we can obtain the 
first order correction as 



" -*P(n').(18) 



>[o], 



E [ ?}(n')-Ef{n) 



The condition with which we can safely neglect the first 
order correction is 



(^ [ °}(n')\0(n',n)\^(n)) 



JO]/ 



E [ l(n')-Ef(n) 



< 1, k' f ,n' ^ k P ,n. (19) 



III. DISSIPATION OF CENTER-OF-MASS 
MOTION 

We consider a moving particle coupled to a boson bath, 
the interaction between the particle and the bath can be 
described by [H[ 



H T =x <g> ^.9fc(6fc + b 



(20) 



where x is the coordinate operator for the particle, 
and b' k are the annihilation and creation operators for the 
bath, respectively. By a standard procedure, we can de- 
rive a master equation and write it in the form of Eq. ([1} , 
but the dissipation term C in this case takes 



£-9 — 7i (2xpx — px 2 
+ 72 (xpp + ppx - 



-x 2 p) 
xpp — ppx) . 



(21) 



V, 



Define C = C T = ^x and D = = ^x - -%L 

the dissipation term can be rewritten as Cp = CpD^ + 
DpC^ - ±{L>tc, p } - i{C*tD,pk which can not be ar- 
ranged in the Lindblad form 5l|. Following the same 
method as that in the last section, we obtain an effective 
Hamiltonian similar to Eq. (JSJ) , but in this section, we re- 
place Ls by Lc with Lc — ^[2-fix 2 — ^i{xp + px)] — 



;[27i(* A ) 



A\2 



„A,„A 



„A„ A)] 



X 



i[2"/ixx 



2L ~ ;iV ^ , "/ 2 (x~P~+P 
,,A 1 We choose the eigenfunctions of Hj 



A 



j 2 (xp A 

px")\. We choose the eigentunctions ot lig = Hs 
as the bases, which are a direct product of the eigenfunc 
tions of Hs and ifg. Let Hs and H A have the eigen- 



functions \xm{x)) and \Xm( x )) wrtn the same eigen- 
values s m for x — x A , we can expand the eigenstate 
of H T by {\A mn (x,x A )) = \Xm(x))\ X A (x A ))} as |$) = 
Smn ^mnlAmn). Following the same derivation pre- 
sented in SecJITl we obtain the following equation for Q mn 
as, 

H T {mn)^ mn + F T (mn)<£ mn + ^ T (mn,pq)$ pq 

pq^mn 

+ L^{mn,pq)<P pq = E$ mn . (22) 

pq^mn 



where 



H T (mn) 



[Vf - iXp(mn)] 



2M 



+ V T (r) + E mn (r) + l£(mn), 
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and 

Lj;(mn,pq)<& pq = (A mn \L c \A pq ). 

When the two subscripts in Lq are same, we will 
omit one of them to shorten the notation, pq ^ mn 
in the summation means (p, m) and (q, n) can not 
be taken equal simultaneously. Other terms are 
similar to the results given in Seed as F T (n) = 

- m 'Epq^mn ( ( A ™« I ^ P A P q ) ( A P? I V>A m n ) ) , T (mn, pq) = 

-2i 7 (2(A mn |V r -A P9 >V r - + (A m „|V|A M )). We would 
like to note that because of the bases we choose are 
uncoupled between the quantum system and ancilla, 
the Hamiltonian H T (mn) governing the center-of- 
mass motion of the system can be factorized as 
flT(mn) = -5m[V - iA{m)f + V +j m - (-^ A - 
iA A (n)] 2 + V A + s n ) + L^(mn) with A(m) = i(x m |Vxm) 
and A A (n) = i(x A \^ A Xn)i where V A represents the 
derivative with respect to x A . This means that the inter- 
action between the system and the magnetic field makes 
the system and the ancilla feel vector potentials A(m) 
and A A (n), respectively. DissCOM induces a correction 
L^ ; (mn) to the system via the interaction between the 
system and ancilla in H T (mn). In the same way, namely 
treat the terms F T , T and L^(mn,pq), (mn ^ pq) as 
perturbations, we obtain the validity condition for the 
zero order BO approximation 

($f, ] (>?)l (0 T (pq,mn) + Ll(pq,mn)) |$L 0] (mn)) 

Ef(pq) -Ef(mn) < ^ 

mn, kf ^ pq, k' f , (23) 

where 

Lc(mn,pq)$ pq = (A mn \L c \A pq ), 

and 

L c = -jp'lix 2 - J2(xp +px)} 



-i[2~/ 1 xx A - j 2 (xp A + px A )]. 

From Eqs. (|22|) and (j23|) we can find that the zero- 
order Hamiltonian H T only contains the diagonal ele- 
ments of Z/£, while the perturbation matrix includes the 
spin part O t , F T , and the spatially dependent term 
L T Q(mn,pq),(mn ^ pq). Recall that the main concept 
of BO approximation in our case is to separate the spin 
variables from the spatial one (i.e., the fast and slow vari- 
ables). Due to the dissipation represented by L^, this 
separation of those variables becomes difficult, resulting 
in the violation of BO approximation. We will show in 
the next section that this is different for the case of spin 
relaxation, because the spin relaxation is coordinate in- 
dependent. 

IV. EXAMPLE 

As an example, in this section we consider a neutron 
moving in a static helical magnetic field, 

- / 2ttz 2ttz \ 

B = B(z) = B (sin 6 cos — , sintfsin — , cos 1(24) 

The Hamiltonian for such a system is 

-2 

H = H(z) = ^- i +pB-S = H K + H s . (25) 

For fixed but arbitrary z, the interaction Hamiltonian 
Hs has the eigenf unctions 

|xi( , ))= ^c 0S |exp(-^)>j 

lx 2 W)=( Sinf ! X c P ( & 7^ } ), (^) 

and corresponding eigenvalues £1,2 = ±//-B. This system 
is the same as that studied in Ref . [22I] . Taken the spin 
relaxation into account (see Eq.([3])), the effective Hamil- 
tonian related to the spin reads 



Hj(^, V A )=pB 



sin(9e~ 



sin f/e' 



2 cos 6» - 




\ l 9 



suit 



sinl 

-2 cos( 
— sin 6 







" ¥9 

,-i<p A 





sin 9e~ lv 
-sm6e lv>A 




(27) 



where <p = 2 jf- and ip A = 2w £ . The rescaled coupling Ej of Hj((j), 4> A ) in this case are given by (in units of /iB) 

constant is defined by g — and therefore it is dimen- 

sionless. In the following, we set 6 = f . The eigenvalue _ 1 

Ei = -7:91, 



5 



£f + \giE) - ~(8 + .g 2 )^ + 2 5 i( e - i ^-^ - 1) - 0, 



= 2,3,4), 



the corresponding right eigenstates are 

Ah 




while the left eigenstates read 



1 



Here 



For j = 1, 



( L j\ = -J" > 



Nj = (Ajaj + Bjbj + C jCj + Djd,). 

Ai = £»i = ai = di = 0, 

B 1 = e-^ A , C 1 = e Ufi , 



bi = e 



up ' 



c\ = e 



-up 



and for j = 2, 3, 4, 



A j =4-giE j -2E% ) 

B 3 = -2e~ i,p gi + 2e^ K E jl 

Cj = 2e lipA gi-2e iv Ej, 



4e ^ A ) +g(g-2iE j ), 



(4 - giEj - 2E]) 



bj = 2e iv E. 

Cj = —2e~ 
d 3 = 4. 



Ej, 



We can see from the eigenvalues and eigenstates that only 
when if — if A , i.e. z = z A , there exists a steady state for 
the quantum system (47j , which is independent of initial 
state. This is exactly the situation we will with in the 
following discussion, in which we will discuss the popula- 
tion transfer among the internal states for the quantum 
system. Initially, the spin of the neutron is prepared in 
the state |+^), we manipulate the neutron moved from 
z = 0toz = Lina time interval T and calculate the po- 
larization of the neutron along the z axis at time T. The 
polarization of the neutron along z axis at time T with 
different rescaled dissipation rate is plotted in Fig[TJ Sev- 
eral features can be found from these figures. (1) When 
g — ► 0, the relation between the polarization along z 
axis and the duration T is a cosine function. This coin- 
cides with the results given in Ref . [22j , meaning that our 
description can return back to the results for closed sys- 
tems. (2) As g increases, the relation between the polar- 
ization along z axis and duration T is also an oscillating 




T(n/\iB) 

FIG. 1: The polarization of the neutron along the z axis as 
functions of the duration T and the rescaled dissipation rate 
g. The time T is plotted in units of tt/^,B. Initially the spin 
is in the state |+|), i.e. P z = 1. 



-g=0.5 — »— g=1 — •— g=2 — — g=20 ^^g=80 




FIG. 2: (Color Online)The polarization of the neutron along 
the 2 axis for different g as a function of gT (in units of tt/^iB). 



function, but the amplitude of the oscillation decreases 
with the interval T. (3) With g -> 0.5 and T -> 3(tt/^S), 
the polarization along z axis reaches zero, arrives at a 
steady value. These features all result from the compe- 
tition between the spin-field coupling represented by H$ 
and the dissipation C. When T and g are large enough, 
the accumulation of the dissipation leads the spin of the 
neutron into the steady state, a result of the balance be- 
tween the spin-field coupling Hs and the dissipation. We 
plot the polarization along z axis as a function of gT for 
different g in Figj2] We can see from the figure that when 
g — * oo, P z sharply drops into -1, i.e. the spin dissipates 
into the state |— i). This can be easily understood as 
follows. When g — >■ oo, the spin-field coupling Hs can be 
neglected, namely the dissipation dominates the dynam- 
ics of the spin, which makes the spin relax to its ground 
state | — i). 

Finally, we discuss the validity condition for the BO 
approximation. For the example under discussion, if we 
do not consider the dissipation, i.e., the system is closed, 
the validity condition given in this paper returns to the 
validity condition derived in Ref.[22|, this means that 
the BO approximation for open systems defined in this 
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0.8 
. 0.6 
0.4 
0.2 



" 



0.2 0.4 0.6 0.8 1 

G=1-e 9 

FIG. 3: Validity measure T{g) as a function of the rescaled 
dissipation rate g. Other parameters in the figure are set to 
satisfy —g^rz = 10 ~ 6 > JSml = 2 x 10-4 • The results have 
been normalized in units of T(g — 0). 



paper can back to closed systems. In order to see the 
effect of dissipation on the validity condition, we recall 
that without any dissipation, the validity condition for 
this example is [221] 



have normalized the results with respect to T(g = 0), i.e. 
all the values in the figure are rescaled by T(0). We can 
find that as g increases, T(g) decreases. This result tells 
us that if the closed system satisfies the BO approxima- 
tion condition, the corresponding open system satisfies 
that condition too. Two points are worth addressing. (1) 
This conclusion depends sharply on the dissipation, for a 
given dissipation, we should check the validation condi- 
tion case by case, and (2) the dissipation benefits the BO 
approximation in this example, this can be understood 
as follows. In this example, the BO approximation for 
open systems is defined as the possibility to separate the 
spin and spatial variables, the dissipation occurs only for 
the spin (slowly changed) variable and it is coordinate 
independent, so it is reasonable that the dissipation ben- 
efits the BO approximation. Similar conclusions can be 
found in the adiabatic approximation in open systems, 
see Ref. [13, El S3. 



CONCLUSION AND DISCUSSION 



1 



fj,BM 2 L 



«1, 



fiBML 



« 1, 



where k z = Namely, when the magnetic field is ho- 
mogeneous and strong enough as well as the velocity of 
the neutron along the z axis is very small, the BO ap- 
proximation is a good approximation. In the following, 
we study how the dissipation affects the validity condi- 
tion for the open system. To show the dependence of 
the validity condition on the dissipation rate, we define 
the following function with max taken over all m and n 
except m = n, 



T(g) = max ■ 



{42}{n)\0{n,m)\<bf{m)) 



■Mi 



E [ °}{n)-E [ ° ] (m) 



(28) 



to characterize the violation of the BO approximation. 
The numerical results of T(g) on g is shown in Fig[3J In 
this figure, we have set all parameters except g to satisfy 
the BO approximation condition for closed system given 
in Ref.H (e.g., = 10" 6 , ^ = 2x HT 4 ), and 



In this paper, we have extended BO approximation 
from closed to open systems by the effective Hamilto- 
nian approach. Two types of dissipation are consid- 
ered, the conditions under which the BO approxima- 
tion holds are given and discussed. When the spin re- 
laxation/dissipation is considered, the center-of-mass of 
the system feels a complex vector potential depending on 
the interaction between the spin and the magnetic field 
as well as the dissipation rate. The zeroth order BO ap- 
proximation was defined as the regime in which its effec- 
tive Hamiltonian can be diagonalized with fixed slowly 
changing variables. When there is only DissCOM, the 
vector potential can be factorized into two parts, which 
are for the system and ancilla, respectively. The Diss- 
COM introduces a correction to the Hamiltonian for the 
center-of-mass motion and then modifies the dynamics of 
the system. An example with spin relaxation/dissipation 
has been presented and discussed. 

We thank Zheng- Yuan Xue for help. This work is sup- 
ported by NSF of China under grant No. 10775023. 
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